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There is a large body of literature investigating the static and dynamic properties of polyelectrolytes due both 
to their widespread application in industrial processes and their ubiquitous presence in biology. Because of 
their highly charged nature, polyelectrolytes tend to alter the local dielectric permittivity of the solution within 
a few nanometers of their backbone. This effect has, however, been almost entirely ignored in both simulations 
and theoretical work. In this article we apply our recently developed electrostatic solver based on Maxwell’s 
equations to examine the effects of the permittivity reduction in the vicinity of the polyelectrolyte. We first 
verify our new approach by calculating and comparing ion distributions around a linear fixed polyelectrolyte 
and find both quantitative and qualitative changes in the ion distribution. Further simulations with an applied 
electric field show that the reduction in the local dielectric constant increases the mobility of the chains by 
approximately ten percent. More importantly, variations in the local dielectric constant lead to qualitatively 
different behavior of the conductivity. 


I. INTRODUCTION 

Electrophoresis is the directed motion of an object in 
an aqueous solution subject to an external electric field. 
Electric fields are often used for the characterization, pu¬ 
rification, and manipulation of polymers^^^, colloids'^^^, 
and cells®’®, which all tend to ionize in aqueous solutions. 
While it is tempting to understand electrophoresis as the 
balance between the external electric force and hydro- 
dynamic friction, the physics is somewhat more compli¬ 
cated^®. The reason for this is the aggregation of oppo¬ 
sitely charges ions, termed counterions, near the surface 
of the object. The ions closest to the surface make up 
the Stern layer, strongly adsorbed ions which co-move 
with the electrophoresing object, effectively reducing the 
charge of the object. Ions further away from the sur¬ 
face form the diffuse layer, which has a characteristic size 
given by the Debye length. Ad: 

/ N 

Ad = JttAb 

V 

where N is the total number of charge species, rij and Zj 
are the number concentration and valency of species j. 
The Bjerrum length Ab is given by: 

Ab = e®/(47reoersfcBE), (2) 

where eg and e^s are the vacuum and the relative permit¬ 
tivity, respectively. 

The ions in the diffuse layer move in the opposite direc¬ 
tion of the analyte since they are oppositely charged. The 
counterions reduce the velocity of the electrophoresing 
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object, through hydrodynamic and electrostatic coupling 
to the analyte. The plane separating these two layers is 
called the shear plane, and the two layers are often called 
electrical double layer (EDL). 

Simulations of polyelectrolyte electrophoresis have, for 
the most part, focused on the electrophoresis of a single 
chain in bulk solution (free solution electrophoresis). The 
first results in this regard were lattice-Boltzmann simula¬ 
tions by Grass and Holm^^“^^, and multi-particle collision 
dynamics simulations by Erank and Winkler^®’^®. Both 
methods were able to reproduce the sharp rise in elec¬ 
trophoretic mobility as a function of the polymer length 
for short chains, followed by a slight decrease for longer 
chains, as is seen in experiments^. In contrast, when 
long-ranged hydrodynamic interactions were ignored, the 
simulations showed a monotonic decrease in the mobil¬ 
ity in stark contrast to experiment. The reason for the 
monotonic decline in mobility in the absence of hydrody¬ 
namic interactions is that a fraction of the counterions 
strongly bond to the polymer backbone, co-moving with 
the chain, and effectively reducing the linear charge den- 
sityS®>s®48_ This phenomenon, termed counterion con¬ 
densation, was first proposed by Manning^®. 

A sharp increase in the mobility with increasing chain 
length is observed for short chains. This increase can 
be understood through the cooperative shearing of the 
fluid within the Debye layer®. Because of this, the mo¬ 
bility starts to flatten out when the chain length becomes 
comparable to the Debye length^®. The aforementioned 
condensation of counterions and subsequent reduction of 
the effective charge of the chain is responsible for the 
slight decrease in mobility for long chains®®’^®’®®. 

A number of studies have also looked at the elec¬ 
trophoresis of polyelectrolytes through simple channels. 
Smiatek and Schmid investigate the role of surface slip 
on the mobility of a polyelectrolyte electrophoresing be¬ 
tween parallel plates®®’®^. More recently, Yan et al looked 
at the electric field dependent mobility of a polyelec- 
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trolyte between parallel plates^^. The electrophoretic 
motion of chains through entropic traps^^ has also been 
studied extensively using computer simulations^^^^^. A 
number of studies have also looked at the electrophoresis 
of DNA in post arrays^’^®’^®. 

Electric helds are also often applied to polyelectrolyte 
solutions in order to measure their conductivity and char¬ 
acterize the properties of the solution^®’®°“®^. Theory is 
able to correctly predict the conductivity under a variety 
of solvent conditions and salt concentrations^®’®®’®^. The 
conductivity normalized by the concentration of charge 
carriers is called the equivalent conductivity and is often 
used instead of the raw conductivity. Deviations from a 
simple linear scaling of the conductivity with respect to 
the concentration, due to both electrostatic and hydro- 
dynamic interactions, causes the variations in the equiv¬ 
alent conductivity. In the salt-free case, there are a num¬ 
ber of experiments which show an initial decrease in the 
equivalent conductivity, followed by an increase start¬ 
ing at polyelectrolyte concentrations of approximately 
c « O.OIM®®^®^. Scaling theories®®’®® are able to cor¬ 
rectly predict the initial decay, but fail to explain the 
increase in the equivalent conductivity at polyelectrolyte 
concentrations above c « O.OIM®®. This increase in the 
equivalent conductivity appears to be independent of the 
type of polyelectrolyte used®®^®^, the length of the poly¬ 
electrolytes®®, and the temperature of the solvent®®. 

In most theories and simulations on charged systems 
the presence of water molecules is taken into account 
as a constant dielectric background. This is most of¬ 
ten done by setting a constant relative permittivity of 
the system Ej. = 78.5, corresponding to salt-free water at 
10°C. However, it has been shown that the presence of 
charged salt ions®®’®® and charged surfaces®® can consid¬ 
erably alter the local dielectric constant. Recent theoreti¬ 
cal work has shown that taking into account the smoothly 
changing dielectric properties between a spherical colloid 
and the surrounding fluid causes a depletion in counteri¬ 
ons near the surface®®^®®. This effectively increases the 
thickness of the Debye layer and in turn reduces the hy¬ 
drodynamic and electrostatic coupling of the surface to 
its counterions. Presumably, a similar effect occurs for 
polyelectrolytes, though the magnitude of the effect is, a 
priori, unclear. Only a few techniques exist for calculat¬ 
ing ion distributions in the presence of a varying dielectric 
constant®®^®®', and for the most part these methods have 
only been used for calculating static ion distributions. 

The influence of a varying dielectric background con¬ 
sists of three separate contributions: The dielectri¬ 
cally reduced pairwise electrostatic interaction between 
charges can be described by the aforementioned algo¬ 
rithms, as well as the influence of dielectric enclosures 
on a point charge by means of virtual boundary or im¬ 
age charges®®’®®^®®, or via a functional approach®®’®'’’®®". 
However, in addition to these two contributions, mobile 
ions have self energy fluctuations in an inhomogeneous 
dielectric medium. In uniform dielectric media, the sol¬ 
vation energy of charged particles is invariant and can 


be absorbed into the chemical potential. However, this 
term cannot be disregarded when the solvent has a vary¬ 
ing dielectric permittivity. As pointed out in an earlier 
publication®®, this direct Coulomb term can be described 
by a modification of the Born energy®® by introducing the 
local self energy contribution. The self energy of ion i is 
given by 


TjseXi _ 

* 47reo 

where the first term including Gpoi is the contribution 
from virtual surface or image charges, and the second 
term relates to the solvation energy within the medium. 
Here, Qi is the charge and Ui is the Born radius of the ion. 
If the gradient of the locally varying permittivity e(r) is 
non-zero, this energy leads to a force that is included in 
the HIM method presented in our earlier work with Xu®® 
and the electrostatic algorithm used in this manuscript. 

In this article, we expand on our previous study on 
the conductivity of polyelectrolyte solutions with a lo¬ 
cally varying permittivity®®. We hrst introduce the new 
electrostatic algorithm in detail. The algorithm extends a 
lattice-based electrostatic algorithm based on Maxwell’s 
equations. We first examine the case of a line of fixed, 
charged monomers and compared the results to Poisson- 
Boltzmann theory around a cylinder. We then apply the 
algorithm to a free polyelectrolyte in aqueous solution 
and examine the influence of the resulting permittivity 
gradients on the couterion distribution and radius of gy¬ 
ration. Next, we look at how the variations in the local 
dielectric constant affect the electrophoretic mobility of 
polyelectrolytes and the conductivity of polyelectrolyte 
solutions. We end with a brief conclusion. 
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II. SIMULATION METHOD 
A. Molecular Dynamics 


We perform Molecular Dynamics (MD) simulations us¬ 
ing the Extensible Simulation Package for Research on 
Soft matter, ESPResSo®®^®®. There is a purely repulsive 
Weeks-Chandler-Anderson (WCA) potential®® between 
all particles to represent steric interactions: 


kwcA(r) 



r < 2®/®aLj, 
otherwise, 

( 4 ) 


where r is the distance between the particles, ctlj = 
0.3 nm is the fundamental MD length scale, and Emd = 
feuT = 4.11 X 10“®® J is the fundamental MD energy 


scale. 

Mp polyelectrolytes consist of a linear chain of nega¬ 
tively charged monomers qi = — le (where e is the fun¬ 
damental MD unit of charge) of length N. Adjacent 
monomers are linked together using a hnitely-extensible 
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nonlinear elastic (FENE) potential 

VfeneW = In ^1 - ^ , (5) 

where Rq = 1.5ctlj is the maximum extension of the 
bond and k = SOeMD/cr^j is the energy scale of the 
bond. An equal number MpN of counterions with charge 
<7 = le are added to the box to keep the net charge of 
the system zero. All particles have a mass m = Imo, 
where nio is the fundamental MD unit of mass. Veloc¬ 
ities and positions are updated using the velocity Ver- 
let algorithm with a time step A^md = O.OIt, where 
T = y/woO’Lj/e is the MD unit of time. In simulations 
where we measure static properties, we use a Langevin 
thermostat with a temperature k^T = £md and a fric¬ 
tion constant FLangevin = jT. All of our simulations 

used periodic boundary conditions in all directions. 


B. Hydrodynamic Interactions 

In simulations where we calculate dynamic properties, 
we chose the lattice-Boltzmann (LB) method to model 
the hydrodynamic interactions, since it is extremely ef¬ 
ficient and has been shown to produce physically sound 
results^’^^’®^’®®. Specifically, we use ESPResSo’s D3Q19 
lattice-Boltzmann (LB)®® with a kinematic viscosity v = 
0.8(t£j/t and density pfluid = Iwo/f^Lj- The LB grid 
spacing is olb = 0.4 nm and the time step was set to 
AfuB = O.Olr. These parameters will give the correct 
hydrodynamic radius of 0.44a « 0.15 nm for counteri¬ 
ons while keeping the lattice-Boltzmann algorithm sta¬ 
ble. We applied an external electric field of Eext = 
O.IeMD/o’Lje (which has been shown to be low enough 
to still be in the linear response regime^®’®^) when mea¬ 
suring mobilities and conductivities. The temperature 
was held constant using an LB thermostat with a tem¬ 
perature k^T = £md- 


C. Electrostatics 

Electrostatic algorithms in molecular dynamics typ¬ 
ically deal with variations in the dielectric con¬ 
stant using a functional, boundary charges, or mirror 
charges®®’^®’'^®’^®^®®’®®’®®. These methods have a number 
of disadvantages for studying polyelectrolytes. Most im¬ 
portantly, such methods do not allow for gradual changes 
in the dielectric constant, but only sharp discontinuities 
in the permittivity. In addition, there is often no way to 
have charges move across a gradient in the dielectric per¬ 
mittivity, or to change the background bulk permittivity 
with time. Our algorithm^® does not have these limi¬ 
tations. It is based on the local electrodynamic solver 
introduced by Maggs and Rossetto in 2002^^. The algo¬ 
rithm was later adapted for off-lattice Molecular Dynam¬ 


ics simulations by Rottler and Maggs^®’^®, and in parallel 
by Pasichnyk and Diinweg^^. 

Traditional approaches typically calculate the electro¬ 
static potential $ by solving the static Poisson equation 

V£V<i> = —p, (6) 

where £ is the dielectric permittivity, and p is the charge 
density. Our algorithm instead focuses on solving the 
electrodynamic equivalent. Gauss’ law 

VD = p, (7) 

where D = eE includes the dielectric permittivity and 
the local electric field E. This works because electro¬ 
statics is merely the limit of electrodynamics where the 
speed of light c approaches infinity. The wave propaga¬ 
tion speed c can be reduced significantly and the algo¬ 
rithm still reproduces the same particle dynamics and 
thermodynamic observables^®^^^. This is very similar 
to Car-Parrinello Molecular Dynamics^®, where the ve¬ 
locity of the electrons is reduced to unrealistic values, 
yet the dynamics of the corresponding nuclei remains 
accurate. In our electrostatic algorithm, the reduction 
of the propagation speed makes the method computa¬ 
tionally feasible and has been shown to require a similar 
amount of computation time compared to more estab¬ 
lished methods^®. The use of electrodynamics for solving 
electrostatic interactions in molecular dynamics simula¬ 
tions has been coined Maxwell Equations Molecular Dy¬ 
namics (MEMD). 

We have already implemented several other algorithms 
that allow for variations in the local dielectric con¬ 
stant within the ESPResSo software package®®’®®’™. The 
method method presented here nevertheless offers a few 
key advantages compared to other electrostatics algo¬ 
rithms. Most importantly, the method is intrinsically 
local, and allows for an arbitrary spatial and temporal 
variation of the local dielectric permittivity. We have 
successfully applied this implementation to a colloid with 
a smooth radial change of the dielectric permittivity in 
the vicinity of the colloid surface®®. We showed that a re¬ 
gion of varying dielectric permittivity, into which charges 
can enter, gives rise to a strong force in the direction of 
the permittivity gradient and therefore significantly in¬ 
fluences the structure of the electric double layer (EDL) 
around the colloid®®. This force is the result of a finite 
sized object in a permittivity gradient, that results from 
the self-energy. This self-energy term is missing from 
methods which only include a discrete jump in the di¬ 
electric constant and significantly affects the equilibrium 
distribution of ions. We also compared our implemen¬ 
tation to a Monte Carlo code that can deal with radi¬ 
ally symmetric and smooth permittivity gradients and 
found excellent agreement®®, giving us confidence on the 
reliability of the algorithm for spatially varying dielec¬ 
tric permittivity. The MEMD grid spacing was set to 
omemd = 0.4 nm, the time step was A^memd = AtMD, 
and the artificial mass was /mass = 0.05mo. 
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III. THE INFLUENCE OF CHARGES ON DIELECTRIC 
PERMITTIVITY 

In MD simulations of polyelectrolytes, an implicit wa¬ 
ter model is often used since explicit water molecules 
would increase the necessary compute time by more than 
an order of magnitude^^. To get reliable results in the 
case of electrophoresis, the implicit water model needs to 
correctly reproduce two properties of water: the hydro- 
dynamic interactions (which we deal with using lattice- 
Boltzmann as explained in section IIB) and the screen¬ 
ing of electrostatic interactions. Most MD simulations 
deal with the screening of electrostatics by introducing a 
bulk dielectric permittivity s. This shows up as a con¬ 
stant prefactor in the Poisson equation (6), and since it 
does not depend on any local parameters, can simply be 
included in the electrostatic force calculation using the 
same prefactor. 

This electrostatic screening effect originates in the re¬ 
arrangement of water molecule dipoles in the presence of 
an electric field E. On average, the dipoles will have a 
tendency to align with the electric field lines and create 
an additional field in the opposite direction. The result¬ 
ing electric field E{r), at a position r, will be weakened 
by a factor of e(T’), which leads to the dehnition and use 
of the displacement field D = e{r)E. As implied by 
the notation e(r), this prefactor is not necessarily con¬ 
stant throughout the system but can depend on the local 
surroundings. Specifically, if there are charged particles 
like polyelectrolyte monomers or salt ions present, the 
water molecules in the vicinity of these charges will be 
restricted in their rotational degrees of freedom by the 
strong electric held create by the charge particle. The wa¬ 
ter molecules thus cannot react to external electrostatic 
inhuences as freely as bulk water molecules. This means 
that the dielectric response in the presence of salt ions is 
reduced. Several studies have examined the dependence 
of the bulk permittivity on the salt concentration in wa- 
^gj.38,78-80^ as well as the variation in the permittivity 
around a charged object^°. 

Most current electrostatics algorithms only allow for 
a change in the system-wide dielectric permittivity, and 
this parameter can be adjusted in polyelectrolyte simula¬ 
tions in an attempt to account for the salt concentration 
in the solution, an approach we also investigate in this 
manuscript. However, it is apparent that the permittiv¬ 
ity is not constant in a system that includes a highly 
charged object surrounded by counterions, and a bulk 
phase that contains mostly solvent and relatively few co- 
and counterions. This is particularly true for a dilute 
polyelectrolyte solution. The local ion concentration will 
not only vary in the proximity of the polyelectrolyte, but 
also in time since the polyelectrolyte is mobile. In the 
present article, we will introduce a method to deal with 
these spatial and temporal changes in the local dielectric 
constant and examine their influence on the static and 
dynamic properties of polyelectrolyte solutions. 


IV. ITERATIVE APPROACH 

As a starting point, we consider a polyelectrolyte to 
be a charged rod, or in our case a linear arrangement of 
charged beads fixed in space as shown in Fig. 1. The rod 
consists of TV = 80 monomers which are fixed in place for 
the duration of the simulation and spaced 0.3 nm apart 
in a cubic box with a side length of 24 nm with periodic 
boundary conditions. Inside the cylinder representing the 
polyelectrolyte, there are very few water molecules and 
the dielectric permittivity is set to £ = Outside 

the polyelectrolyte, we expect the dielectric permittivity 
to depend on the local concentration of counterions, and 
gradually increases to around e = 78.5 in the bulk. 



FIG. 1. Schematic of our first simulation setup: a fixed line of 
monomers (blue) is placed in the center to form a stiff rod and 
counterions (red) can move freely around the polyelectrolyte. 
The permittivity (yellow background) is set to e = 2 within 
the rod and is variable outside the rod, eventually reaching 
e = 78.5 (white) at longer distances. 

If we average the salt concentration over time, we 
should get a cylindrical ion profile around the rod cen¬ 
ter because of the system’s symmetry. This amounts to 
a one-dimensional problem to solve, and there are cer¬ 
tainly faster and more efficient ways to approach it than 
an MD simulation. However, we wanted to introduce a 
verification for our algorithm, and we wanted a point of 
reference for later comparison when the dielectric con¬ 
stant varies both spatially and temporally. To this end, 
we will solve for the electrostatic forces using the imple¬ 
mentation of MEMD with varying dielectric permittivity 
in ESPResSo. 

To calculate the local permittivity value from the lo¬ 
cal ion concentration, we require a function that maps a 
given salt concentration to a permittivity value. Eor our 
simulations, we use the empirical function obtained by 
Hess et al.®® from atomistic MD simulations of a sodium 
chloride salt solution: 

78.5 

^ “ 1 -h 0.278 • C ’ ^ ^ 

where C is the salt concentration in moles per liter [M]. 

We self-consistently solve for the counterion distribu¬ 
tion and the dielectric constant using an iterative scheme. 
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We start out with a flat permittivity profile, which as¬ 
sumes a constant e = 78.5 for the surrounding solvent. 
We then map this salt concentration to a dielectric per¬ 
mittivity via equation (8), and interpolate the resulting 
permittivity distribution to a grid. We then set the per¬ 
mittivity values of the lattice links and, after equilibra¬ 
tion, run the simulation until the average radial counte¬ 
rion distribution converged. We use a successive under¬ 
relaxation, where the two preceding ion distributions are 
averaged to stabilize the iterative scheme. The averaged 
ion distribution is then used to generate new permittivity 
values and the process is repeated iteratively until both 
the ion distribution and permittivity converge. The re¬ 
sults of the first eight runs are presented in Fig. 2. 



0 2 4 6 8 10 
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FIG. 2. Starting from a flat permittivity profile of e = 78.5, 
the counterion distribution including a salt concentration de¬ 
pendent permittivity around a charged stiff rod is calculated 
iteratively. The counterion distribution converges to a sta¬ 
ble solution (Run 8, top graph) which is significantly differ¬ 
ent from the often applied Poisson-Boltzmann solution (thin 
black line, top graph). The first iterative run expectedly gives 
a result identical to the Poisson-Boltzmann distribution. The 
dielectric permittivity, calculated from equation (8) including 
an additional contribution from the polyelectrolyte charge, 
shows a weaker attraction between the counterions and the 
rod (bottom graph). 

The simulations clearly converge towards a stable final 
distribution, and for the simulation shown in Fig. 2 it is 
converged after the eighth iteration. The converged coun¬ 
terion density actually increases in the first two nanome¬ 
ters, before an exponential decrease from two nanometers 
to around five nanometers from the rod. This roughly 
corresponds to the region of reduced dielectric permit¬ 
tivity shown in the bottom graph of Fig. 2, which also 
extends approximately five nanometers from the rod. It 


is counterintuitive at first glance that a decreased permit¬ 
tivity close to the polyelectrolyte would lead to a decrease 
in counterions at the surface, since the electrostatic at¬ 
traction increases. However, the effect seen here is mostly 
due to the force that the permittivity gradient Ve exerts 
on a counterion, which pushes the charges away from the 
backbone, along the permittivity gradient. 

The counterion distribution in Fig. 2 looks similar to 
our findings for spherical colloids^^, which is to be ex¬ 
pected given that the underlying physics is independent 
of the specific geometry. Because of these similarities and 
the very stable final distribution, we are confident that 
our simulation method produces a realistic counterion 
distribution for a stiff rod model including varying per¬ 
mittivity effects. However, this method is time consum¬ 
ing and not suitable for dynamic simulations since the 
dielectric constant does not vary temporally. Nonethe¬ 
less, it does provide a good point of comparison for the 
adaptive approach introduced in the next section. 

V. ADAPTIVE APPROACH 

To gain a more flexible and widely applicable algo¬ 
rithm, we opted to calculate the instantaneous local salt 
concentration during the simulation. We again simulate 
a stiff rod consisting of N = 80 monomers fixed in space 
with a distance 0.3 nm between monomers, in a cubic 
box with a side length of 24 nm with periodic boundary 
conditions. 



FIG. 3. To calculate a local salt concentration on the fly, we 
use a weighted average in the surrounding 7^ lattice cells. The 
weights are equal to the inverse of the shell number, meaning 
the central cell gets a weight of 1, the second shell a weight 
of 1/2^, the third shell a weight of 1/3^ and the outermost 
cells getting a weight of 1/4^. This calculation is carried out 
for every MEMD cell individually. The result can then be 
mapped to a dielectric permittivity for the center cell using 
equation (8). 

For every lattice cell we calculate a weighted average of 
the ion concentration within a predefined cube, as shown 
in Fig. 3. The ion concentration is calculated for ev¬ 
ery single MEMD cell in the simulation box and treats 
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charged monomers exactly the same as the salt ions. This 
means that the local ion concentration is simply a metric 
for how many charged particles are in the vicinity of the 
MEMD cell. If the mesh is chosen too coarsely or the 
screening length too small, there will be large jumps in 
the permittivity when ions enter or leave the cube. This 
should be avoided since dielectric jumps are physically 
unrealistic and represent barriers that block ion flow al¬ 
most entirely. For our system, we have chosen a cube size 
of 7^ lattice cells. Given the fine mesh of the simulation, 
this corresponds to a reach of 1.4 nm or 21b (Bjerrum 
lengths) in water, which is a realistic screening length 
and we found it to result in sufficiently smooth permit¬ 
tivity curves. The charges included in the surrounding 
lattice cells are weighted by the inverse square of the shell 
number (see Fig. 3). This represents a 1/r^ influence, in 
accordance with the assumption that the polarization re¬ 
sponse is linear to the electric field, which also decays as 
1/r^. 

For verification of this new adaptive scheme to calcu¬ 
late the local charge concentration, we again simulated 
the system sketched in Fig. 1. We did not manually set 
the dielectric permittivity within the polyelectrolyte and 
ran the simulation starting with a random distribution 
of counterions. The result and comparison to the iter¬ 
ative approach for the counterion distribution and the 
permittivity is shown in Fig. 4. 



FIG. 4. Comparison of the newly adopted scheme to the final 
distribution reached with our iterative approach. The inset 
shows the permittivity distribution around the rod. Both 
results are almost identical. 


The solution provided by our new adaptive scheme re¬ 
sults in almost the exact same structure of the electric 
double layer. Instead of a monotonic exponential de¬ 
crease, the counterions are pushed away from the poly¬ 
electrolyte and reach a stable maximum at around 1.1 nm 
from the surface. This shows that the counterion distri¬ 
bution is not sensitive to the exact method used to cal¬ 
culate the local dielectric constant, which indicates that 
our somewhat ad hoc method of calculating the local ion 
concentration correctly captures the underlying physics. 


VI. FLEXIBLE POLYELECTROLYTE 

With this new adaptive scheme, we are now able to 
fully simulate a flexible polyelectrolyte in aqueous solu¬ 
tion with the local dielectric constant dependent on the 
local salt concentration. The simulation is still compu¬ 
tationally demanding, since the MEMD lattice has to be 
very fine to accurately represent the changes in dielec¬ 
tric permittivity around the polyelectrolyte, and because 
of the calculation of the salt concentration for every cell 
using the surrounding 7^ neighbor lattice sites. It is, how¬ 
ever, far less costly than simulations with explicit water 
molecules. Depending on the number of processors, it 
runs about a factor of 2 to 3.5 slower than simulation 
with constant background permittivity. 

For a flexible polymer, it is not as straight forward 
to obtain a radial distribution function of counterions. 
However, for better comparison we adopted the calcula¬ 
tion presented in Fig. 5. 



FIG. 5. To calculate the radial distribution function of 
counterions for a flexible polyelectrolyte, we form a cylinder 
around each monomer, the axis of which is parallel to the 
connection between the two direct neighbor monomers. 

For each monomer, we calculate the vector between 
the centers of the two neighboring monomers to which it 
is bonded. This vector is then shifted to the monomer’s 
center. Finally, we calculate the radial distribution func¬ 
tion in a cylinder of length 2ci, where d is the equilibrium 
distance between two monomers. 

It should be pointed out that with this scheme, it is 
possible for a single counterion to be included in the dis¬ 
tribution function of several different monomers at the 
same time step. To ensure a valid distribution function, 
it is therefore imperative to keep a running total of the 
total number of counterions counted throughout the sim¬ 
ulation and normalize the distribution function accord¬ 
ingly. 

To verify that the scheme works with a flexible 
polyelectrolyte, we simulated a single chain with 50 
monomers, each carrying a charge of (jm = ~le. To keep 
the system neutral we added 50 counterions with charge 
Qci = le but no additional salt ions. We again use a pe¬ 
riodic box with side length 24 nm. A comparison of the 
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result for this flexible polyelectrolyte to the schemes ap¬ 
plied to the stiff rod, to a simple Poisson-Boltzmann solu¬ 
tion, and to atomistic simulations are presented in Fig. 6. 
The three models that include a salt dependent variation 
of the dielectric permittivity show quite good agreement. 
They all differ qualitatively from the Poisson-Boltzmann 
solution. Note that we previously showed that simula¬ 
tions with a sharp dielectric jump turns out to be almost 
identical to the Poisson-Boltzmann result®®. 



(a) comparison of the three methods 



(r-R) [nm] 

(b) atomistic simulations 

FIG. 6. Comparison of the counterion distribution functions 
for a flexible polyelectrolyte (red triangles), a stiff rod us¬ 
ing the iterative (dashed line) and adaptive (dotted line) 
schemes, and the Poisson-Boltzmann solution that does not 
take smoothly varying permittivity into account (solid line). 
The depletion of counterions close to the backbone is present 
in atomistic simulations as well. The initial rise in the counte¬ 
rion concentration and thickening of the Debye layer is due to 
a force pushing the connterions away from the polyelectrolyte 
backbone. This force is proportional to the gradient of the di¬ 
electric constant and physically is the result of the solvation 
energy of an ion being higher in a dielectric background with 
a larger dielectric constant. 

An additional confirmation is that our simulations 
qualitatively match atomistic simulations, as seen in 
Fig. 6b. The all-atom Molecular Dynamics simulations 
were performed with the GROMACS 4.5.5 software pack¬ 
age®® at 300 K for a model polyelectrolyte in aqueous 
solution and in the presence of counterions. The sys¬ 
tem studied is a sodium chloride (NaCl) solution with 
a Kirkwood-Buff based force field®^ in combination with 


the SPC/E water model®®. Kirkwood-Buff force fields 
have been shown to reproduce thermodynamic and static 
properties in good agreement to experimental findings 
and to avoid the spurious artifacts of other force fields, 
like the overestimation of ion pairing effects. We con¬ 
structed a simple fictitious linear polyelectrolyte with 30 
’CH2’ beads as defined in the GROMOS force field®® 
with the corresponding monomeric binding distance of 
0.149nm. The polyelectrolyte can be interpreted as a 
rod of infinite length by using the periodicity of the sim¬ 
ulation box with a cubic side length of 3.804 73 nm in 
agreement to the approach presented in Ref.®^. We as¬ 
signed a charge of q = ±le to every second monomer 
while the other monomers remain uncharged, which gives 
a line charge density of I = ±3.94 e/nm. The Bjerrum 
length Ab = 0.78 nm for the SPG/E water model®^ yields 
a Manning parameter C, = IXb/c = 3.08 ^ I which in¬ 
dicates a large fraction of condensed counterions. We 
randomly inserted Nc = 15 counterions in the box to 
achieve electroneutrality, which gives a salt concentration 
of c = 0.45 M. Finally, the positions of the monomers 
were fixed which avoids the influence of configuration ef¬ 
fects on the ion distribution®®. The largest difference 
between the coarse-grained and atomistic simulation in 
Fig. 6 is that the counterion density decays much faster in 
the case of the atomistic simulations. This is because the 
density of counterions in the box is much higher, which 
effectively decreases the Debye length Ad®®- 

A. Radius of Gyration and Diffusion of an Isolated 
Polyelectrolyte 

Despite the relatively large difference in the counterion 
distribution around the chain when variations in the local 
dielectric constant are taken into account, most static 
properties of the polyelectrolyte, like the radii of gyration 
or end-to-end distance, were surprisingly unaffected by 
including the local drop of the dielectric constant in the 
vicinity of the polyelectrolyte. This is clearly seen in the 
radius of gyration values plotted in Fig. 7. Note that 
here we scaled the box length with the cube root of the 
chain length L = 27A^®/®, in order to keep the counterion 
density constant, and thus the Debye length. Both with 
and without taking into account variation in the local 
permittivity constant, the radius of gyration scales as 
jyo.®, good agreement with previous work®^’®®’®®. 

The diffusion of the polyelectrolyte, as shown in Fig. 8, 
leads to slightly different parameters when fitting the 
power law 

D = Dox-"^ (9) 

to the data, where the fit parameter Dq represents the 
diffusion coefficient of a single isolated monomer. The 
exponent decreases from m = 1.03(±0.02) for constant 
permittivity to m = 0.93(±0.03) for varying permittivity. 
This is somewhat surprising since based on Zimm’s the¬ 
ory we would expect that the diffusion coefficient would 










FIG. 7. The radius of gyration as a function of the chain 
length with a uniform permittivity (black diamonds) and with 
varying permittivity (red triangles). We have included expo¬ 
nential fit functions that give a scaling behavior of 0.79 and 
0.78 for simulations with and without varying permittivity, 
respectively. 


FIG. 9. Normalized electrophoretic mobility of the polyelec¬ 
trolyte for constant and varying dielectric background. While 
both data sets show a similar qualitative behavior, the mo¬ 
bility is significantly increased when variations in the local 
permittivity are taken into account. For comparison, and to 
judge the influence of ghost mode thermalization, we have 
included the results obtained by Grass et al. 


scale with the same exponent as the radius of gyration. 
The fact that the diffusion coefficient drops off slightly 
faster is partly due to the rather small size of the peri¬ 
odic box, which creates large finite-size corrections. The 
main contribution is, however, that the scaling law is only 
strictly valid in the limit of long chains. Fig. 8 shows that 
the scaling exponent 0.8, obtained from the radius of gy¬ 
ration data, can also reasonably be used to fit the data 
at large values of N. 
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thermalized ghost modes and finally to a new simulation 
featuring locally and temporally varying permittivity, as 
calculated by our adaptive approach. Our simulations 
were carried out using a single polymer of lengths N, with 
one charge per monomer. N counterions with a 

charge qd = le were also added to the system, but no ad¬ 
ditional salt. The box size was set to L = 27iV^/^(TLj to 
maintain a constant counterion density (which also keeps 
the effective Debye length constant®®) for all simulations. 
We measured the velocity of the polymers in the direction 
of the field and divided by the magnitude of the external 
field to get the mobility values in Fig. 9. The mobility is 
normalized using /ijed = /^/mo with /rg = e/QinjlB- The 
inclusion of spatially and temporally varying permittivity 
shows the same qualitative behavior, but produces mo¬ 
bilities which are significantly larger than the simulations 
with a uniform dielectric constant. 

The mobility values can be interpreted as the ratio of 
an effective charge to that of an effective friction coefh- 

cient3.ii-i6.i8 


FIG. 8. Normalized diffusion of a flexible polyelectrolyte with 
and without varying dielectric permittivity. Both simulations 
follow the predicted D = Dox~^ behavior. The scaling pa¬ 
rameter changes slightly from m = 0.93(±0.03) for varying to 
m = 1.03(±0.02) for constant permittivity. However, this is 
mainly due to deviations from the expected m = 0.8 behavior 
at very short polymer lengths, as can be seen by comparing 
with the dotted line. 


B. Electrophoresis 

In this section, we compare the data of Grass and 
Holm^^ including hydrodynamic interactions to the same 
simulation with a Lattice-Boltzmann algorithm including 


= 


Qeff 


( 10 ) 


where Qeff and Fes are the effective charge and effec¬ 
tive friction of the chain respectively. We can calculate 
the effective charge per monomer of the polyelectrolyte 
three different ways: (i) Dynamically, by measuring the 
mobility in a Langevin dynamics simulation in response 
to an externally applied electric held, (ii) Statically, by 
determining the integrated charge shift of the moving 
inhection point according to the ion distribution. And 
(iii) we can calculate a theoretical prediction with the 
assumption of Manning condensation. 

According to Manning’s theory, the predicted value for 
condensed counterions for a Poisson-Boltzmann distribu- 
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tion is 


cThondlB ^ 0.243 nm 

“ « 0.714 nm 


1-0.34, (11) 


where crLj.bond is the average bond length, and Ib is the 
Bjerrum length. Manning’s theory is known to provide 
good predictions of the effective charge in the limit of 
long chains, but fails for very short chains^^’^®’^®. 

Regardless of the chain length one can calculate the 
(static) effective charge directly using the method devel¬ 
oped by Belloni®^ and Deserno®^. We plot the integrated 
charge density as a function of the logarithm of the dis¬ 
tance to the backbone in Fig. 10. The counterions closer 
to the backbone than the inflection point of this graph 
are considered condensed^^’®^. 


second derivative 



log(r-R) 

FIG. 10. Normalized integrated charge density around the 
polyelectrolyte, plotted against the logarithm of the distance. 
The inflection point of the graph, pointed out by gray dashed 
lines, can be determined via a second derivative (upper plot, 
the dotted line marks zero) and is seen as the boundary be¬ 
tween condensed and free counterions. For varying permit¬ 
tivity, this inflection point moves (green arrow) further away 
from the surface and to a lower integrated charge value, sug¬ 
gesting a higher effective charge and mobility of the polyelec¬ 
trolyte. 


In Fig. 10 we see that including variations in the local 
permittivity causes the inflection point to move (green 
arrow) to longer distances r from the backbone. At the 
same time, the number of condensed counterions Ned 
decreases, increasing the effective charge QeS = {N — 
Ned), resulting in an increase electrophoretic mobility of 
the polyelectrolyte. 

In the case of Langevin dynamics simulations, one can 
calculate the effective friction Fes directly, resulting in: 


Qeii 


Qei^ 

ivf ’ 


( 12 ) 


where N is the number of monomers, and F = ImojT 
is the friction coefficient of an individual MD bead. A 
rescaled version of the mobility, Qes/N = F/r, is plotted 
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FIG. 11. Effective charge of the polyelectrolyte calculated in 
three different ways. The dashed line in both graphs depicts 
the effective charge predicted by Manning^®, (a) shows the 
dynamic effective charge calculated via the measured mobil¬ 
ity with an applied electric field and without hydrodynamic 
interaction (Langevin dynamics), (b) shows the static effec¬ 
tive charge calculated from the integrated charge density as 
sketched in Fig. 10. The simulations with (red triangles) and 
without (black diamonds) varying dielectric permittivity show 
very little deviation (gray squares, right axis) at shorter poly¬ 
mer lengths and the difference increases until N ~ 15. 


in Fig. 11a and matches Manning’s prediction in the long 
chain limit. 

Both the results for the Langevin simulations and the 
static effective charge calculations in Fig. 11 yield very 
similar results, with a difference of less than 3% in the free 
draining limit. For a constant dielectric background, the 
effective charge for polyelectrolytes longer than Af = 15 
monomers fits well with Manning’s prediction^®. Our 
simulations with a flexible varying dielectric permittivity 
show a significant change compared to our simulations 
with constant permittivity in effective charge, with the 
difference being roughly 15% for long chains. This ex¬ 
plains why we see an approximately 15% increase in the 
electrophoretic mobility when taking into account varia¬ 
tions in the local dielectric constant. 

We should add that taking into account only an in¬ 
crease in ionic attraction due the lower dielectric constant 
in the vicinity of the polyelectrolyte, as done for exam- 
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pie in Ref.®^, typically leads to a stronger electrostatic 
attraction, and therefore to a decrease in the effective 
charge of the polyelectrolyte. However, in our simula¬ 
tions the dielectric permittivity changes gradually in the 
vicinity of the polyelectrolyte, which results in a solva¬ 
tion force that pushes counterions further away from the 
backbone. This leads to a significantly widened electric 
double layer, as shown previously in Ref. and in sec¬ 
tion IV, and therefore to less tightly bound counterions 
and an increased effective charge. 

We have shown that the inclusion of varying dielectric 
permittivity leads to a significant increase in the elec¬ 
trophoretic mobility as well as the number of free counte¬ 
rions, since fewer counterions are condensed to the poly¬ 
mer backbone. These two effects should combine in an 
additive way when looking at the conductivity. In our 
simulations including an external electric field, we mea¬ 
sured the conductivity cr using the relation; 


where J is the current density, and E is the applied ex¬ 
ternal electric field. 



N 


6? 


T) 


FIG. 12. The conductivity as a function of the polymer length 
N. The deviation between simulations with and without vary¬ 
ing permittivity is even more pronounced, since the higher 
mobility is accompanied by more uncondensed charges that 
can contribute to the conductivity. 


The conductivities with and without variation in the 
permittivity were converted to SI units using the relation 






(gMD)2(^SI)3^MD’ 


(14) 


where ctlj is the length scale, e is the unit charge and C is 
the friction coefficient of an ion. The superscript denotes 
the according units system. To calculate the friction co¬ 
efficient of an ion in SI units we used the Stokes relation 
= 67r77(T®\ where 77 = 8.9 x 10“"^ Pas is the dynamic 
viscosity of water. For the friction coefficient we 
used the expression from Ahlrichs and Diinweg®^ 


c 


MD 


1 

f 


1 1 
lA 2577a’ 


(15) 


where Tq (the bare friction), 77 (the dynamic viscosity), 
and a (the grid spacing) are the lattice-Boltzmann pa¬ 
rameters. The results are shown in Fig. 12 and are very 
comparable for short polymer lengths, but for > 12 
taking into account variations in the local dielectric con¬ 
stant results in a 20% increase in the conductivity. To 
relate this to experimentally accessible results, we are in¬ 
terested in the ratio between the measured conductivity 
and the ideal conductivity The ideal conductiv¬ 
ity simply assumes that the mobility of all charge carriers 
is /i = q/67r77i?Hj which essentially ignores the hydrody¬ 
namic and electrostatic coupling between charge carriers. 
This lets us define a correlation coefficient 


A = 1 - 


.M 


rid ’ 


(16) 


which is 0 when there is neither hydrodynamic nor elec¬ 
trostatic interactions between the charge carriers. The 
correlation coefficient quantihes the amount of friction 
and coupling between the polyelectrolyte and the coun¬ 
terions. The ideal conductivity can be determined with 
the diffusion constants of both species, via the Stokes- 
Einstein equation 


Ne^D+ , N'^e'^D- 
"" ksT + ksT 


(17) 


where we can measure the diffusion constant for a single 
particle and determine the diffusion constant of the 
polyelectrolyte with 


1 

D- = -^- 

N GTrrjRH 


(18) 


The hydrodynamic radius Rh in this equation is defined 
as 


1 

Rh 



(19) 


which we calculated directly from our simulations. The 
results for the correlation coefficient A are plotted in 
Fig. 13. 

We observe a significant difference between the ap¬ 
proach with and without dielectric variations, exceed¬ 
ing 30% for longer polyelectrolytes. This clearly shows 
that the coupling between the polyelectrolyte chain and 
the surrounding counterions is significantly reduced when 
variations in the permittivity are taken into account. The 
difference is the result of the structural differences in the 
EDL. This setup should be more readily accessible ex¬ 
perimentally, and allow for a quantitative comparison to 
our simulations. 


C. Conductivity for varying monomer concentration 

We now turn to the conductivity of polyelectrolytes 
in salt-free aqueous solution. In Fig. 12, the conduc¬ 
tivity plateaus at polymer lengths above 20 monomers. 
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FIG. 13. The correlation coefficient shows a significant de¬ 
crease for the simulation with varying dielectric permittivity. 
This is due to the weaker coupling between the polyelectrolyte 
and the counterions, and should be observable in experiment. 


With the same simulation setup as before, we measured 
the conductivity for different monomer concentrations by 
keeping the box size constant at (32 nm)^ and adding 
polyelectrolyte chains and counterions. We calculated a 
mean value and standard deviation for the conductivity 
with lengths N = 30, 45, and 60 to make sure that the re¬ 
sults are independent of N. Renormalizing this data with 
the monomer concentration gives an equivalent conduc¬ 
tivity per monomer A. We also normalize all the data to 
an extrapolated value of A(C' = 0) = 1. Through this 
rescaling A is independent of both the friction coefficient 
of the MD beads and the viscosity of the fluid. 

Let us first consider the case of constant background 
permittivity, the blue line in Fig. 14. Note that we set 
the background permittivity using Eq. 8, although the re¬ 
sults are almost identical to simply using a constant per¬ 
mittivity of e = 78.5 for all monomer concentrations®®. 
The equivalent conductivity drops dramatically until a 
monomer concentration C « 0.01 M. In this regime, the 
Debye layer shrinks leading to more and more counteri¬ 
ons condensing on the polyelectrolyte backbone as seen in 
Fig. 15. At higher monomer concentrations C > 0.01 M, 
both the equivalent conductivity and fraction of con¬ 
densed counterions continue to increase, but at a much 
slower rate. 

Simulations including variations in the local permit¬ 
tivity also display an initial decrease of the conductivity 
with increasing polyelectrolyte concentrations, which is 
again due to an increase in the fraction of condensed 
counterions in Fig. 15. There is a clear minimum around 
C = 0.01 M, after which the equivalent conductivity be¬ 
gins to rise. This can be attributed to the decrease in 
the fraction of condensed counterions in Fig. 15. The 
decrease in the fraction of condensed counterions can in 
turn be related to the large decrease in the local relative 
permittivity from around 55 in the dilute limit, to ap¬ 
proximately 41 at the highest monomer concentrations 
(the red line in Fig. 15). This must be the result of 
an increase in the local ion concentration, since it alone 
determines the permittivity in our simulations via equa- 



Colby, Na D Kwak, Na ^ constant e(C) 

Kwak, Cs O Kwak, Li ^ varying s 

Kwak, K ® Lipar-Ostir 



C [M] 


FIG. 14. The equivalent conductivities A/Aq as a function of 
the molar monomer concentration C. The experimentally ob¬ 
served minimum (gray symbols) is reproduced in simulations 
with varying e{r) (red line), while simulations with constant 
but scaled e{C) (blue line) exhibit a different qualitative be¬ 
havior. 
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FIG. 15. Average permittivity in the vicinity of the poly¬ 
electrolytes £poiy as calculated by the adaptive scheme (black 
squares, left axis), and the fraction of condensed counterions 
fed (right axis) for constant but scaled e{C) (blue circles) and 
varying e{r) (red triangles), fed has a maximum for variable 
permittivity around the conductivity minimum observed be¬ 
fore. The salt dependent permittivity keeps decreasing even 
though the number of condensed counterions almost stagnates 
at higher concentrations, which is due to the coiling of the 
polyelectrolytes. 


tion 8. The reason for the increased ion concentration is 
that the decreasing Debye length causes the polymers to 
coil significantly as can be seen in figures 16 and 17. 

We also compare our results in Fig. 12 to existing ex¬ 
perimental data from Kwak and Hayes®'^,Colby et al. 
and Lipar-Ostir et al. Since the raw results greatly 
depend on the hydrodynamic radius of the solvated coun¬ 
terions and differ by more than a factor of 2, we normalize 
all the data to an extrapolated value of A((7 = 0) = 1. 
The experimental data show excellent agreement with 
our simulation results. The slight differences are most 
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(a)C = 0.001M (b)C = 0.01M (c)C = 0.1M 


FIG. 16. Snapshots from our simulations. (a) At low 
monomer concentration, the polyelectrolyte has few con¬ 
densed counterions (red) and is stretched out because of elec¬ 
trostatic repulsion of the backbone, (b) At higher monomer 
concentration, the number of condensed counterions is signif¬ 
icantly increased, (c) At high monomer concentration, the 
polyelectrolyte coils. In the simulations with varying dielec¬ 
tric permittivity depending on the local charge concentration, 
this coiling of highly charged monomers leads to a lower av¬ 
erage permittivity around the polymer backbone. 



FIG. 17. Average radius of gyration Rq of the polyelectrolytes 
as a function of the polyelectrolyte concentration for a uni¬ 
form background permittivity e{C) (red squares, solid line) 
and a locally varying permittivity e{r) (black circles, dashed 
line). 


likely due to enthalpic factors between the counterions 
and the polyelectrolytes or the specific structure of the 
polyelectrolyte, however, these are clearly secondary fac¬ 
tors. 


VII. CONCLUSIONS 

We used our novel implementation of a local electro¬ 
static method to include spatial and temporal changes 
in the local dielectric permittivity. We then performed 
simulations of polyelectrolytes solutions. We first looked 


at a fixed charged straight polymer, and found excellent 
agreement of the counterion distribution with an addi¬ 
tional iterative method. Overall, the results were quite 
similar to findings for the counterion distribution around 
colloidal particles^^’"'^. We have also showed that the 
inclusion of a dielectric jump in the dielectric permittiv¬ 
ity at the surface of the polyelectrolyte was insufficient 
to even qualitatively reproduce the same counterion dis¬ 
tribution®®. This is because the dielectric jump fails to 
take into account the force caused by the gradient in 
the dielectric permittivity related to the solvation energy 
of the ions. Our results show that this force can not 
be neglected in this setup, since it significantly widens 
the electric double layer and therefore reduces counte¬ 
rion condensation. While many studies have looked at 
the effect of reduced permittivity within a polymer, col¬ 
loid, or surface, our results clearly show that the gradient 
in the dielectric constant adjacent to the interface actu¬ 
ally plays a much larger role. 

We then applied our algorithm to a free, flexible poly¬ 
electrolyte. Interestingly, the distribution of counterions 
around the backbone of the fluctuating polymer are al¬ 
most identical to the results for the fixed, straight poly¬ 
mer. Despite a qualitatively different structure of the 
electric double layer around the charged polymer, we ob¬ 
served almost no change in the radius of gyration of the 
chains when taking into account salt-dependent varia¬ 
tions in the local dielectric constant. This demonstrates 
that the change in the local counterion distribution only 
has a minor influence on the conformational properties 
of the polyelectrolyte. 

Next, we looked at the role of varying permittivity on 
the mobility of a polyelectrolyte in aqueous solution sub¬ 
ject to an external electric field. We found that the qual¬ 
itative behavior of earlier simulations^^’^® is reproduced, 
in agreement with experimental data^®’^^. That being 
said, the difference in the distribution of counterions 
within the electric double layer had a significant quan¬ 
titative influence on the electrophoretic mobility. The 
mobility and effective charge of the polyelectrolytes, in 
comparison to a simulation without varying permittivity, 
shows an average increase of around 9% for long poly¬ 
electrolytes. We observed a roughly 20% increase in the 
electric conductivity when taking into account changes in 
the local permittivity, since both the increased effective 
charge of the polyelectrolyte and the increased number of 
uncondensed counterions increase by approximately 10%. 

We then investigated the role of monomer concentra¬ 
tion on the conductivity of salt-free polyelectrolyte solu¬ 
tions. There was a clear qualitative difference between 
the simulations with constant dielectric background and 
the ones using our adaptive scheme for variations in the 
permittivity, where the latter does not show a continuous 
decrease of conductivity with rising monomer concentra¬ 
tion, but exhibits a distinct minimum and subsequent 
rise in the equivalent conductivity. The same nonmono¬ 
tonic behavior has previously been observed in experi- 
ments®®^®^, and our simulation results closely match ex- 
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periment and even exhibit the conductivity minimum at 
the same monomer concentration. 

Our results can be succinctly summarized by saying 
that the varying dielectric permittivity causes a thicken¬ 
ing of the Debye layer, which has only a minor influence 
on the conformation of the chain. Despite not making 
dramatic changes to the static properties of the system, 
the electrophoretic mobility and electric conductivity is 
significantly larger. Experimental data for the equiva¬ 
lent conductivity at increasing monomer concentration 
is reproduced using our scheme for varying permittiv¬ 
ity, while the simulations assuming a constant dielectric 
background show qualitatively different behavior. This 
indicates that correctly accounting for variations in the 
local permittivity due to the local salt concentration is 
necessary to achieve quantitative and qualitative agree¬ 
ment with experiment. Most likely taking into account 
local variations in the dielectric permittivity would lead 
to similar quantitative differences for other systems such 
as the electrophoretic mobility of colloids, and electroos- 
motic flow in microfluidics. 
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